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ABSTRACT 


Accurate prediction of human response to sonic booms from proposed HSCT aircraft depends 
on a knowledge of the waveshape and risetime of the boom at the ground. In previous work, we have 
developed a numerical technique to predict the combined effects of molecular absorption and finite 
wave distortion on the sonic boom as it propagates from the aircraft to the top of the turbulent boundary 
layer. We have more recently developed a scattering center based model to calculate the effects of tur- 
bulence on the sonic boom waveform as it propagates through this boundary layer. Calculations have 
been performed using single scales of turbulence and compared to measurements at Edwards AFB in 
the late 1960s. A model of the atmosphere involving two scales each for convective and mechanical 
turbulence has been developed and fit to meteorological data collected during JAPE 2. Scattering calcu- 
lations employing this model undeipredict the number of unperturbed waveforms. In order to develop 
a more realistic model of the atmosphere, the JAPE 2 meteorological data has been fit to a von Karman 
spectrum. Results of scattering using this multi-scale model will be presented. The combination of 
finite wave effects with turbulent scattering predictions includes the principal effects of the atmosphere 
on the sonic boom from the HSCT. 
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INTRODUCTION 


The prediction of the average environmental impact of the HSCT requires accurate modeling of 
the processes affecting the sonic boom waveform and risetime. We have used the enhanced Anderson 
algorithm to predict the risetime and waveshape of sonic booms under non-turbulent conditions. This 
method can also be used to predict the risetime and waveshape at the top of the turbulent planetary 
boundary layer. 

The enhanced Anderson algorithm includes all finite wave effects and the vibrational relaxation 
effects of N2, O2, and CO2 in combination with atmospheric H2O. This algorithm has been compared 
to data from explosions 1 and sonic booms 2 and has been tested against measurements of high intensity 
ballistic waves from rifles and from tank guns 3 . In addition, the results of this calculation for quasi 
steady shocks agree with the results from the enhanced Burgers’ Equation 4,5 . 

Figure 1 presents the results of the application of the enhanced Anderson algorithm to a 
predicted HSCT waveform 6 . We emphasize that the key parameter in determining the risetime of the 
sonic boom is the absolute humidity. 

Under turbulent conditions, the risetimes of sonic booms are scattered and are occasionally as 
large as ten times the risetimes calculated from vibrational relaxation considerations. It is clear that 
turbulence is the cause of the increased risetime and peculiar waveforms observed. Analytic techniques 
have been used to estimate the increase in average risetimes 7, 8 > 9 and to calculate perturbed waveforms 
due to focusing and defocusing of the waves by turbulence 10 . In such calculations, it is usually 
necessary to assume a single strength and turbulence scale representative of the atmospheric turbulence. 
The largest turbulence effects are usually identified when the largest scales are chosen as typical. 

We have chosen a different approach to calculating the effects of turbulence on sonic boom 
risetimes and waveforms based on a simple scattering center-based theory. The scattering center-based 
method accurately predicted the effects of turbulence on the coherence of continuous wave signals 
above natural ground surfaces 11 . 


METHOD 

The scattering center-based technique resolves atmospheric fluctuations into a sum of discrete 
spherically symmetric Gaussian “turbules”. The total effect of the atmosphere is then calculated by 
summing up the scattering amplitudes. See Figure 2 . The scattered amplitudes are calculated using the 
first Born approximation. If the complex pressure at the receiver is written as: 
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where the superscript B refers to the first Bom approximation, p 0 (r) is the unperturbed spherical wave, 
and N is the number of turbules, then 
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s defines the 1/e 2 contour of the turbule, qi is the index of refraction profile strength, and 0 O is the 
scattering angle. The geometry is indicated in Figure 3. 

The initial research on continuous wave propagation modeled the atmosphere as a random sum 
of identical turbules. This single scale calculation was extended to impulse propagation with promising 
results. 12 The impulse is Fourier transformed into the frequency domain and the total scattered compo- 
nent at each frequency is calculated. Then the inverse Fourier transform yields the time domain wave- 
form. The single scale calculation (s = 10m, 30m or 100m) with a fluctuating index of refraction of 
<p. 2 > = 10 x 10' 6 predicted spiked and rounded waveforms and predicted risetimes as large as 10 ms. 

These results encouraged us to analyze the results of the JAPE-2 tests ^44 using the scattering center- 
based model. 


ANALYSIS OF JAPE-2 DATA 

The JAPE-2 tests consisted of simultaneous measurement of sonic boom characteristics and 
meteorological measurements. The wind and temperature fluctuations were measured at heights up to 
30m using sonic anemometers and hot wire anemometers. The sonic boom data was analyzed by 
Wiltshire, Garber and DeVilbiss 1 ^ and provided as computer files. The turbulence data was analyzed 
by Bass, Boulanger, Olsen and Chintawongvanich 15 . 



a.) Two Scale Model 


Examination of the data showed that a single scale model of the atmosphere could not fully 
describe the turbulence above the ground. The time correlation of the fluctuation quantities was fit to a 
two scale model. See Figure 4. Table I displays the results of the analysis for a moderately turbulent 
day during JAPE-2. 

Table I. Example of the Two Scale Model Applied to Atmospheric Data 



Wind driven 1 

Wind driven 2 

Temp, driven 1 

Temp, driven 2 

Size (meters) 

117 

11 

74 

8 

Number of 
Eddies 

1 

90 

1 

233 

(n 2 ) 

0.54 10- 5 

0.25 lO- 5 

0.5 10-« 

0.4 10- 6 


The scattering calculation was performed by summing the results of four calculations - one for 
each scale size. The input waveform to the scattering calculation was an N-wave propagated from the 
flight altitude to the top of the turbulent layer using the enhanced Anderson algorithm. The results of 
the Anderson algorithm agree moderately with the measurements taken under low turbulence condi- 
tions. See Table II. 

Table II. Comparison of Measured and Predicted Waveform Parameters for the T-38 



Measurements for the low 
turbulence case 

Calculations using the 
Anderson algorithm 

Peak overpressure (psf) 

0.71 

0.88 

Risetime (ms) 

0.32 

0.33 


Figure 5 compares the results of the measurement and prediction for T-38 overflights under 
moderate turbulence conditions. Although the scattering center model produces a wide distribution of 
risetimes, it does not predict the shift of the histogram maximum to 2 ms; rather the maximum remains 
at the unperturbed value of 0.3 ms. It is believed that this is due to the use of two relatively large scales 
to represent the atmospheric turbulence. The scattering from large turbules is predominantly in the 
forward direction, and large turbules are relatively sparse, so that it is easy to “miss” the receiver with 
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the scattered wave. The four scale model does, however, represent a significant improvement over the 
single scale model. 

b.) von Karman Spectrum Model 


The fit of the autocorrelation to two scales rather than one improved the prediction of risetimes 
significantly. The high occurrence of unperturbed risetimes indicated that smaller and intermediate 
scales were needed to fully describe the scattering of sonic booms by turbulence. 

De Wolf 16 presented a technique for simulation of a turbulent atmosphere obeying the von 
Karman spectrum in terms of the number density of turbules. 

The general form of a 3-D von Karman spectrum is given in terms of frequency by: 

<t>(f) = af2 ,, (5) 

fc+i# 

where: 

/ \ 2/3 
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The coefficients a and b are determined by fitting a function <()(f) through the measured spectra. See 
Figure 6. 

The fit parameters are then used to determine n (s), the number density of turbules of size s 
needed to model the fluctuating atmosphere. De Wolf’s model was originally developed to predict 
second moments of a scattered field and therefore is designed to reproduce only second moments of the 
fluctuation fields. Higher moments must be accurately represented to express the temporal characteris- 
tics of an impulse. De Wolf used an index of refraction maximum for each turbule of ± 1.0 and em- 
ployed a very sparse distribution. We have varied the product of qi 2 and n (s) until the model distribu- 
tion approximates the measured second and fourth moments <ja 2 > = 9.6 x 10" 6 , <q 4 > = 2.5 x 10* 10 . 
The variation of calculated <p 2 > with number of turbules and qi 2 is shown in Table III. 
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Table IIL Calculate <p 2 > and <|LX 4 > as a Function of Number of Turbules 


Number of 
Turbules 

Percentage of 
Volume 

qi 2 

msm 

W 

42000 

8 

1.5 10- 4 

1.2 10- 5 

1.0 10-9 

63000 

12 

hhbekqbhh 




18 


1.0 10-5 

6.9 10 - 1 ° 

127000 

r 24 

5.0 10- 5 


5.5 10 - 1 ° 

254000 

48 

BMElIEiflHi 

■HEEIESIH 

4.4 10 - 1 ° 


The turbule spatial and size distribution for each realization is determined by Monte Carlo 
methods. The index of refraction fluctuations along a straight line has been compared to the corre- 
sponding measured values and exhibits similar fluctuation scales and displacement. 

The second improvement to the scheme was the use of the measured height of the Planetary 
Boundary Layer in the calculation. Figure 7 displays the temperature versus height curve for one flight 
during JAPE-2. One sounding is taken with the tethesonde going up and the inversion height is 400m, 
the other trace is the tethesonde coming down 30 minutes later and the inversion height is at 670m. The 
turbulent layer thickness at the time of the later sonic boom measurement was extrapolated from this as 
750m. 

The results of this calculation for 20 realizations are displayed in Figure 8. The maximum 
occurrence risetime shows a shift away from the non-turbulent risetime of 0.3 ms. The smaller and 
intermediate scales of turbulence have a significant effect on the risetimes of sonic booms. It is clear, 
however, that the shift is not large enough to match the measured data in Figure 5a. 

CONCLUSION 

The enhanced Anderson algorithm provides a good prediction of waveshape and risetime of the 
HSCT at the top of the Planetary Boundary Layer. 

The scattering center-based model can be extended to predict distorted wave shapes and longer 
risetimes. At this stage, the scattering based model does not predict long enough average risetimes, but 
does show that smaller and intermediate scales are important in increasing the average risetimes. 

The larger scales are the source of the dramatically distorted waveforms, but are not the source 
of the shift in average risetimes. The scattering center-based calculation allows the quantitative investi- 
gation and modeling of the turbulence effects discussed qualitatively by Crow, Plotkin and George, and 
Pierce. 
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Figure 1. a.) Predicted HSCT waveform close to the aircraft. 



b.) HSCT waveform at the ground predicted using the enhanced Anderson algorithm. 
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Figure 5. a.) Measured risetime distribution for the T-38 aircraft. 



Rise Time (ms) 

b.) Predicted distribution using the two scale model. 
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Figure 6. Fit of the von Karman spectrum to a 23 minute sample of wind data. 



Figure 7. Temperature versus height for one flight at JAPE-2 


135 






